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The problem of constructing the kinetic equation with the description of motion of a hole in 
systems with strong spin- hole interaction (such as high- temperature superconductors) in terms 
of the spin polaron has been considered in the framework of the regular antiferromagnetic s — d 
model. It has been shown by the example of the electrical resistivity that kinetics is determined 
by the properties of the bands of the spin polaron (rather than "bar hole") and their quasiparticle 
residues Zk. The cases of low and optimal doping of the Cu02 plane have been considered. It has 
been shown that the rearrangement of the spectrum of the lower polaron band, as well as the strong 
doping dependence of the quasiparticle residues is decisive in the unified consideration of these 
cases. 
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It is known that the normal state of high- temperature 
superconductors is characterized by a complex behavior 
of the spectral and transport properties due to the strong 
interaction of the carriers with the spin subsystem [l|, 0, 
Q . This refers to the nontrivial evolution of the hole and 
spin subsystems with increasing doping, when the system 
passes from the Mott dielectric to the metallic state 

The overwhelming majority of the studies Q on the 
microscopic description of the kinetics in high- temper- 
ature superconductors was devoted to the case of the 
optimal doping and was based on the concept of the 
almost antiferromagnetic Fermi liquid described by the 
spin- fermion Hamiltonian H on the square lattice 

Ht=H + Hf, Hf^-d'E^, 
H = Hoo + J, Hf)Q — Hoh + I, (1) 

Hoh — £kQko.aktT, 

k,(T 

I=lni-p)J2^^+sSn + llpT.^R+'i^^^ (2) 

R,g R,d 

^^7/^ ^ °^k+q,7i'^q'^7i72°'k72 = X^^q'^q' 

^ k,q,7i,72 q 

ak+q,7l'^7l72ak72 (3) 

^ k,7i,72 

The term H^h in the Hamiltonian ifgo describes the 
bar Fermi carriers and contains the spectrum of bar 
holes; / corresponds to the frustrated antiferromagnetic 
interaction between 5=1/2 spins, g and d are the vec- 
tors of the first and second neighbors, respectively; and 
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/i = (1 — p) and I2 — pi, where p (0 < p < 1) is the 
frustration parameter, are the respective antiferromag- 
netic exchange constants. The term J describes the of 
the the carriers with the subsystem of localized spins Sr 
((?" are the Pauli matrices and the summation over re- 
peated Cartesian superscripts a and spin subscripts 71 
and 72 is implied). The total Hamiltonian Ht includes 
the interaction Hf with the electric field E, and is the 
dipole momentum operator. 

In order to adequately describe the temperature de- 
pendence of the electrical resistivity p(T) as well as the 
Hall coefficient Rh{T)) by solving the kinetic equation 
in the almost antiferromagnetic liquid model, the spec- 
trum Sk in Hamiltonian Hoh is always changed to the 
spectrum E^^ of the lowest quasiparticle band of the 
spin polaron, and the operators ato- in Hoh are remained 
fermion operators. The spectrum E^^ corresponds to a 
large Fermi surface and is well measured in experiments 
on angle resolved photoemission spectroscopy (ARPES). 
Such a change is not obvious and inevitably leads to the 
incorrect number of holes rih ~ 1.2 instead of the real 
relation < 0.2. However, for the physical quantities 
appearing in the kinetic equation and in the expression 
for the average current, this change seems reasonable, 
because the kinetic must be determined by the veloc- 
ities of quasiparticles, V^^-* = (?i5^^^/9k, of the lower 
polaron band rather than by the velocities of bar holes, 
Vk = dsi^/dk. 

In this study, we present the derivation of the kinetic 
equation that is based on the Hamiltonian U]) and spin- 
polaron concept and leads to the observed p(T) depen- 
dence for the large Fermi surface and rth < 0.2 (due to 

the small quasiparticle residues Z^^ < 0.2 of the lower 
polaron band in the Green's function of the bar hole), 
as well as to the natural appearance of the quasiparticle 
velocities V^^^ = dE^'>/dk. 
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The carriers are described in the multipole approxi- 
mation, i.e., in the sufficiently complete basis of the spin 
polaron. Under the assumption that doping stimulates 
the frustration of the spin subsystem, the cases of small 
and optimal doping of the Cu02 plane are discussed. 
The simultaneous consideration of these cases is possible 
only with allowance for the strong rearrangement of the 
residue function Z^"' with doping. 

Let us first consider the equilibrium Green's functions 
for bar holes and introduce spin polaron states by an 
example of the simplest approximation. 

At characteristic values J ~ 0.1 eV, the Hamiltonian 
J corresponds to the strong interaction of bar hole with 
the spin subsystem. As a result, elementary charge exci- 
tations must be described by the spin polaron, which is 
represented as a superposition of the bar- hole operator 
Okcr and the spin polaron operators describing "dress- 
ing" of the Okcr to the operators of the spin subsystem. 
The problem is solved with the use of the Mori- Zwanzig 
projection method for Green's functions. The method 
implies the choice of a finite set of the basis operators, 
which must include the pairing of the bar hole with lo- 
calized spins from the very beginning. 

It is known [6] that the minimum " good" site set is the 
following set of the basis operators: 



2^ e ^r+p°^crcri "^r '^ai 



(4) 



(5) 



p,qSf2 

f7 = {q: \±{-Klg)^q^,y\<L}. 



The first two operators ^p^J and Lp\fJ can be treated as 
local spin- polaron operators, the following two operators 
(/3rCT and (/3rCT correspond to the spin polaron of the inter- 
mediate radius and describe the pairing of local polaron 
operators tp^J and tp^^J with the spin wave operators 



(2) 



E 



giq(r+p)_g<Q 



A feature of operators (O is that they reflect the pair- 
ing of spin waves S'q with momenta q close to the an- 
tiferromagnetic vector Q = (tt, tt), i.e., q values fill re- 
gion $7 consisting of four L x L squares in the cor- 
ners of the first Brillouin zone (in what follow, we set 
il ^ L X L = 0.25(71/5)^). Pairing with such spin waves 
S'q takes into account the sharp peak of the spin- spin 
structure factor in the region close to the antiferromag- 
netic vector Q and is responsible for the splitting of the 
lower quasiparticle band appearing in the local polaron 
approximation Q. Moreover, the inclusion of the finite 
region $7 is necessary for describing the correct passage 
to the limit T ^ 0. 



The standard projection procedure for solving the 
equations for Green's functions in the momentum rep- 
resentation for operators (g]) and ^ provides four bands 
of the spin polaron E^^ (s is the band number), an ex- 
plicit expression for the Green's function of bar hole, 
G'h(k, a;) =<C akalajccr ^uj, the expression for the num- 
ber of bar holes, nko- in terms of quasiparticle residues 
Z^P , and makes it possible to represent Hamiltonian H 
given by Eq.([T]) in the polaron basis as Hp: 



Zi 
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nk. = (aL«k.>=E4'^"^(4'^); (6) 



Here 7ij;-(i?k) = (e^^'' + ^, where ^ is the chem- 
ical potential. 



is)' 



(7) 



s s—1 



k,s,7 



(8) 



Here, P is the projection operator on the polaron 
space. The matrix U^^(k) is expressed in the explicit 
form in terms of the spin correlation functions, which are 
in turn determined in terms of the susceptibility x(q, 
of the frustrated antiferromagnetic spin subsystem. 

The Hamiltonian Hp contains the operator Jp, which 
is a component of J that is responsible for the formation 
of polaron from a bar hole. Operator Jp includes only 

those matrix elements of J which describe the polaron 
scattering processes without change in the quasimomen- 



tum, i.e., JpCt^^P —> ct'^^' processes. This scattering 
processes are taken into account in the projection- ap- 
proach description of the formation of polarons. 

Comparison of the projection- method results for the 
a complex spin polaron [6] with the self- consistent Born 
approximation (SCBA) calculations (at T = 0) Q cer- 
tainly indicates that the lower band E^'^ and residues 
ZP''^ well reproduce the SCBA peak and its intensity. 



The upper three bands E^ , E^ , and E^' effectively 
describe the incoherent part of Aincoh{k,u!) of the total 
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hole SCBA spectral function 

AsCBAiK U) = - 4'^) + A,nc^hiK ^) 

The motion of the bar hole under the action of, e.g., 
the external electric field in terms of polaron operators 

(s) 

is the motion of spin polaron a\^^ simultaneously in four 
bands with the velocity vj^*'' = dE^^ /dh. 

The operator J such that H — Hp + J should be intro- 
duced in the kinetic equation. Here, J should be treated 
as an operator including the matrix elements of J, which 

lead to polaron- polaron scattering Ja^^^^ —> a^^^q cr'^-q 
with change in the quasimomentum q 7^ and with the 
simultaneous excitation of the spin subsystem. This ver- 
tex obviously must describe the collision term in the ki- 
netic equation for polarons in the second order in J. 

The Hamiltonian in the polaron representation takes 
the form (from now on, we take term / out of the polaron 
Hamiltonian Hp in explicit way) 



Here, 



where 



H — Hq + J, Hq — Hqo + Jp — Hp + I 



k,i7,S 



(9) 



A" = 
q 



E (E "k+qsl4+q!7lKl72(E 



k.71,72 sl = l 



Ht=H + Hf, Hf = -d^'E, 



(10) 



(11) 



k.k' 



where 



a;kk' = (k|x|k') = (Ak',k+q + A 



(<Z,0); X. 



k.k' 



k',k 



lim(ij=>0)7 

(12) 



k'k 



(s) . . . . 

The transition to operators a^^. is implied in the op- 
erator (P and the quasi- inhomogeneous field E directed 
along the X axis is introduced according to the Eq. pTjl . 

In order to obtain an expression for the electrical re- 
sistivity, we use a variant of the linear- response theory 
in which the field E that ensures a fixed electrical cur- 
rent j (rather than the current at a fixed field) is sought. 
Deviation from equilibrium (k is the system deviation 
parameter specifying j) at the initial time f = is char- 
acterized by a finite set of operators Ff and the density 



matrix p is specified in the form 

+ p^ = Zo'cxp{-Ho/T), 
(A) = Sp{p'^A}, 

l.s 



(13) 



S CIS 



{Fno^O, F, 



I \t=o - 2^ 

k,(T 



^ll(.nFtl), (14) 



[Ff,i7o]=0, [/^r.p^l^O. (15) 
Note that the moments -Ff (k) in the current state are 



1 (i)t (1) 
k"k,- tti--- corre- 



odd in k and the moment F\ — ^ V 
sponds to the simplest one- moment approximation as- 
sociated with the lower polaron band. 

The evolution equation for the density matrix pt — 

]i+ p't has the form = [{Hq + J + Hf), (p + Pt)] 
and its solution for pt is sought with an accuracy of the 
first order in Hf ^ E ^ kA^, assuming that (fP ^ k 

and J ^ X, where A is the scattering parameter. In this 
approximation, the transition to the interaction repre- 
sentation [A{t) = giHot^^-iHot-^ provides 



Pt 



[Hf{T-t),p"]dT+ / [J{T~t),^dT 



+ {-^)^ dr dr'[J(r-i),[J(r'-t),p]]. (16) 



The conditions of the quasistationarity of the current 
density matric p in the limit of infinite time reduce to 
Fflt^aa — Ff\t=Q, i.e., to the system of equations 

Sp{p'tFn\t~.oo = 0. (17) 

In the limit of the infinite number of moments, system 
(|17p is equivalent to the exact kinetic equation for the 
nonequilibrium one- particle density matrix. The equa- 
tions of system (jl7|) have the usual kinetic form 



^ f dTSp{F{[Hf{T~t),p'] 
Jo 



dT / dT'Sp{Fi'[J{T ^t),[JiT' -t) 

a Jo 



2ni 



= 0, 
\l8) 



where the first and second terms correspond to the field 
and collision terms, respectively. We denote these terms 
as tXi and tPf, respectively. Equations (fTT]) determine 
coefficients 77* in (jfi. 

The detailed form of the collision term 
in Eg. psp includes expression of the form 
Sp{FfA-^ iT)Sl (f )A^, iT')SP^ (r')p°Fr/ }> where 
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T = T — t and t' = T — t, which are calculated 
with the use of the mode- coupling approximation 
In this approximation, "outer" operators 

[i.e., operators entering into J in Eqs. PU)) ] are 
separately averaged at the first stage and the initial 
averaging is performed for them. The unaveraged 
averages of the remaining operators that appear af- 
ter the above procedure are calculated at the next 
stage: 5p {F,^A^^ (f (f )A^^ (f')^^, (r')p°F// } => 

[Sp FfA^^ (?) A^, (r')/5°Fr/ ) {Sp p'S^ {r)S^, (?')) 
As a result, for the collision term, we obtain 

pr=2^5:c^;r' (19) 

PuT-J"^ E (^^Ak)-Ff(k + q)) 

^nF{E\:\l-nF{E\:i\)] 

xnB{E\:l^ - 4" V"(q, - 4''^)) (20) 

where x"(Qi ^) is the imaginary part of spin susceptibility 
and nB(a;) is the Bose function. 

Under the assumption that in the expression for (fp 
in Eq. is quasidiagonal, because Hf in Eqs. pTjl is 
quasidiagonal, the field term Xi in Eq. is modified 
to the form 

xf = i^5:y,^(^)fz^^^')z(^)pr(k). (21) 

Finally, the average current is obtained in the form 
X ^ r;fFf (k)l/^^(^)z^n^(4^^)[l - np{E\:\ (22) 

l,sM 

where is the distance between Cu02 planes (we take 
az ~ 6.6 A and the corresponding volume of the unit cell 
is g^a^ ^ 93 A^). 

With the use of the coefficients ?7f obtained from equa- 
tion Xf + = (fT8|) . the current density and diagonal 
component of the resistivity tensor p can be determined. 

For the doping case of interest, rih < 0.2 the chemi- 
cal potential /i lies sufficiently deep in the lower polaron 
band. For this reason, only the lower band s = 1 can be 
retained in the expressions for the current and field and 
collision term; thus, the summation over s is removed in 
Eqs. ([20land ((22l Thus, the problem reduces to a usual 
one- band case, where the spectrum of the lower polaron 
band provides the characteristics of the carrier spectrum. 
The significant difference from traditional expressions is 



that the residue function Z^^ explicitly enters into ki- 
netic equations (^0)1 - ([^ and into Eq. ^ for Uh. Only 
the lower polaron band will be taken into account and its 
index 1 will be omitted in the notation of energy, residues, 
velocities, and coefficients t] and moments f;(k) in Eqs. 

It is known that the dynamics of charge carriers in 
Cu02 planes is well described by the three- band Emery 
model [1, . The calculation of the spin polaron spec- 
trum with the use of the Emery model [6*] provides the 
spectrum observed in the ARPES experiments in a wide 
doping range. The assumption of the correspondence be- 
tween doping in models with free carriers and frustra- 
tion p in the pure spin model [ll| is a key assumption 
for the description of the properties of the lower polaron 
band. This assumption is physically natural: a moving 
hole destroys the magnetic order and the same effect oc- 
curs with increasing p in the pure spin model. Moreover, 
it is based on the similar character of changes in spin cor- 
relation functions as functions of x and p. However, strict 
statements on the x p correspondence are absent, but 
frustrations is always present in the spin subsystem at 
the doped Cu02 plane. Even in the dielectric limit, the 
ratio of exchange on the second neighbors to exchange 
on on the first neighbors is estimated as /2//1 ~ 0.1 
(12j. The role of frustration as a driving force of the for- 
mation of various spin- liquid states is widely discussed. 
It is expected that a quantum phase transition can oc- 
cur near /2//1 ~ 0.5 (which corresponds to p w 0.3) fl|. 
Close frustration parameter values are accepted when dis- 
cussing the stripe scenario of the appearance of incom- 
mensurate peaks Q- 

A decrease in the spin correlation length C, with in- 
creasing p corresponds to change in the spin correlation 
functions, which explicitly appear in the equations for the 
Green's functions of the spin polaron and significantly af- 
fect the spectrum _Ek and residue function Z\^. Below, we 
will discuss the problem of the electrical resistivity for the 
cases A and B of a small spin correlation length on the 
order of several lattice constants C, :^ g {p ~ 0.3) and a 
large spin correlation length Q > lOg {p ~ 0.1), respec- 
tively. Cases A and B refer to cuprates close to optimal 
doping and to strongly undoped cuprates. 

Spin fermion Hamiltonian ([1]) applied to the Emery 
model for case A provides the spectrum E^ and residue 
function Q, whose characteristic form is shown in 
Figs. 1 and 2, respectively. The analytical form of E^ 
is approximated with the use of the square- symmetry 
harmonics 

=^(ai7<,(k) + a27g(k) 
+a37d(k) + a47d(k) + a57g(k)7d(k)); (23) 

where 

7g(k) = (cosgfc^ + cosgky)/2, (24) 
7rf(k) = cos gk^ cos gky, (25) 

and fli — 1.5, 02 = 3.0, as = —1.25, 04 — 0.0, and 
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FIG. 1: Fig. 1. Hole spectrum (1231) presented by contours 
= const in units t = 0.2 eV in the first quarter of the 
Brillouin zone. 




FIG. 3: Fig. 3. Residue function (low doping) presented 
by contour curves in the first quarter of the Brillouin zone. 
The thick line is the Fermi surface in case B. 




FIG. 2: Fig. 2. Residue function Z^^ (doping is close to 
optimal) presented by contour curves in the first quarter of 
the Brillouin zone. The thick line is the Fermi surface in case 
A. 



as — 0.1. The following energy- parameter values char- 
acteristic of the Emery model are assumed below: r = 0.2 
eV, J = T, I = 0.5t, and temperature T = 0.1/ « 120 
K. 

The form of is shown in Fig. 3. The form of the 
spectrum qualitatively corresponds to the following 
scenario of the evolution of -Ek with increasing (: the 
bottom of the spectrum (see Fig. 1) is shifted along 
the r ^ M diagonal toward the point r(0, 0). This shift 
at low doping ensures a "large" Fermi surface located 
near the magnetic Brillouin zone X X. This scenario 
is confirmed both by ARPES experiments and by calcu- 
lations of the spin polaron spectrum. In this case, the 
spectrum E^ near the Fermi surface and magnetic Bril- 
louin zone is inevitably more flat than the spectrum E^. 
This flattering can be treated as the narrowing of the 
band with increasing the polaron effect (decreasing frus- 



tration). In order to taken into account this narrowing, 
we set E^ = 0.75 E";^. The electrical resistivity is calcu- 
lated for the Fermi surfaces shown by thick lines in Figs. 
2 and 3 for case A {fi^ = 0.82t) and B (^^ = 0.375r), 
respectively. 

In order to explicitly clarify the role of the flattering 
of Ek and the k dependence of Zk, we use the following 
spin susceptibility for both cases 

X(q,c^)= , . ■ (26) 



This form of the spin susceptibility appears in the frame- 
work of the spherically symmetric self- consistent ap- 
proach [l^, [3| in the method of irreducible two- time 
retarded Green's functions {{S^\S?^q)}^+ie [H, [3 , [11 
or in the memory function method [l7|. The spin sus- 
ceptibility is taken in the same form as in the case of 
the calculation of the electrical resistivity at T = 0.1/ in 
where the method of self- consistent calculation of 
x(q, oj) and choice of damping 7 = 0.45/ were discussed 
in detail. 

To take into account strong scattering anisotropy due 
to the strong scattering of carriers from the spin mode 
with the antiferromagnetic vector Q, it is inevitably nec- 
essary to go beyond the framework of the traditional one 
moment approximation when solving the kinetic equa- 
tion. The large number of moments is also necessary for 
demonstration of the convergence of the method. 

Polynomials of the velocity components and their 
derivatives will be used below as the moments //(k) of 
the distribution function: 



dkv ' 



dkv 



dk'^ dky 



k ; 



V^}. (27) 



For good convergence, it is sufficient to take into ac- 
count first three or four moments. 
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Calculations for case A give the electrical resistivity 
Pa — 85.2 pVLcm and the number of holes — 0.21. This 
value is comparable with p{T — 120i^) w 100 p,Vtcm for 
La2-xSr^CuOi with 0.16 < a; < 0.22 4]. The electrical 
resistivity pA calculated under the assumption = 1 is 
denoted as 'pA (below, all the quantities calculated under 
the assumption = 1 are marked by tilde) and is = 
SA.2pD.cm ~ Pa and « 1.13. The values pA and 
Pa are approximately equal to each other, because Z^ 
depends only slightly on k in the k- space Brillouin- zone 
region that it is near the Fermi surface and contributes 
to kinetics. In this region, Zk ~ const = Z . In this case, 
r]i ~ Z~-^ according Eqs. (fT9|) - (|2T|) and the resistivity 
p determined by the current (f22|) is independent of Z. 
In this case, the results obtained in the spin- fermion 
models, where the residue function is disregarded, are 
valid. However, the smallness of Z^ {Z^ < 0.2), as well 
as the k dependence of in the entire k- space Brillouin 
zone, see Fig. 2, obviously leads to the relation < 
0.2 at a large Fermi surface and affects the position of 
the Fermi surface with respect to the magnetic Brillouin 
zone. The last effect significantly determines the collision 
integral for strong scattering by the antiferromagnetic 
vector Q. 

In case B, and the Fermi surface shown in Fig. 3 
correctly represent the properties of the spectrum of un- 
derdoped cuprates: the hole residues Z^ decreases from 
0.24 to 0.04 when moving along the Fermi surface from 
the point of intersection with the — M to thre point of in- 
tersection of the Fermi surface with the boundary X — M 
of the Brillouin. This decrease qualitatively reflects the 
known opening of the pseudogap on the Fermi surface. 



Calculations for the case B lead to pB = 231.1 pricm and 
the number of holes rih = 0.08, which is close to p{T — 
120K) « 220 pncm for La2-xSrxCuOi with x « 0.1 
With Zk = 1, we obtain pB — 162.1pncm fa pA and 
Ufi ~ 0.92. Thus, the inclusion of the k dependence of Zk 
becomes important in the case of low doping. Moreover, 
the empirical approximation ri/i sa 1 — n/i is completely 
violated. 

Our calculations also demonstrate the strong depen- 
dence of p on the band narrowing, which is described 
by the Drude formula p — m* /ne^r for simple metals. 
Indeed, the electrical resistivity calculated with the un- 
flattened spectrum (the rigid band, but with residues 
Z^, see Fig. 3) for the Fermi surface position in the Bril- 
louin zone shown in Fig. 3 is equal to p'b = 124.9 pflcm 
(<C 231.1 /ificm); in this case, remains unchanged: 
nil ~ 0.08. The calculation for the case B (we recall that 
the Fermi surface in case B is always located as in Fig. 
3) and the rigid band i^k = in the "rigid" residue 
approximation Zk = Z^ gives lower electrical resistivity 
p'b = 89.5 pVLcm. The calculation for the rigid band and 
Zk = 1 provides p'g = 82.2.5 pflcm and ~ 0.92. 

Thus, we not only have derived the kinetic equation 
for the spin- polaron carriers, but also have shown that 
both change in the residue function Zk and band nar- 
rowing must be included in the description of kinetics on 
the basis of the Fermi surface obtained from the ARPES 
measurements for various doping degrees (e.g., the Fermi 
surface shown in Figs. 2 and 3). 
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